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Measuring central-spin interaction with a spin bath by pulsed ENDOR: 
Towards suppression of spin diffusion decoherence 
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We present pulsed electron-nuclear double resonance (ENDOR) experiments which enable us to 
characterize the coupling between bismuth donor spin qubits in Si and the surrounding spin bath 
of ^^Si impurities which provides the dominant decoherence mechanism (nuclear spin diffusion) at 
low temperatures (< 16 K). Decoupling from the spin bath is predicted and cluster correlation 
expansion simulations show near-complete suppression of spin diffusion, at optimal working points. 
The suppression takes the form of sharply peaked divergences of the spin diffusion coherence time, in 
contrast with previously identified broader regions of insensitivity to classical fluctuations. ENDOR 
data suggest that anisotropic contributions are comparatively weak, so the form of the divergences 
is largely independent of crystal orientation. 
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I. INTRODUCTION 

Quantum decoherence presents a fundamental limitation 
to the realization of practical quantum computing and of 
other technological devices which actively exploit quan- 
tum phenomena. In 2002, a ground-breaking study estab- 
lished the usefulness of so-called optimal working points 
(OWPs):^ These are parameter regimes where the system 
becomes - to first order - insensitive to fluctuations of exter- 
nal classical fields. We consider here the effect of OWPs in 
a system where decoherence of a central spin system arises 
from interactions with a fluctuating bath of surrounding 
spins - a scenario that is of considerable significance in the 
field of quantum information^^ — 

A promising approach for silicon-based quantum- 
information processing (QIP) involves combined electron 
and nuclear spins of donor atoms in Si, which are amenable 
to high-fidelity manipulation by means of electron param- 
agnetic resonance (EPR) and nuclear magnetic resonance 
(NMR), respectively. Most studies have considered phos- 
phorus (^^P) donors in SiJ^— More recently, several dif- 
ferent groups have investigated another Group V donor, 
209 g- These studies not only showed that bismuth donors 
have similar properties to Si:P, such as long electron spin 
coherence times T2 of the order of several ms at low temper- 
atures (< 16 K)^i2i but that they also offer new possibili- 
ties for QIP. For example, strong optical hyperpolarization 
was demonstrated ^i2^ allowing for efficient initialization 
of the nuclear spin. The SiiBi spin system has an elec- 
tron spin S* = 1/2 and a large nuclear spin / = 9/2 as 
well as an atypically strong hyperfine coupling constant, 
A/2t: = 1.4754 GHz. The strong state-mixing occurring 
for magnetic fields 5 c^ 0.1 — 0.6 T where the hyperfine 
interaction competes with the electronic Zeeman energy 
allows transitions which are forbidden at high magnetic 



fields^i21 observed recently in Ref. [H. In Refs. ^^ a 
set of minima and maxima were found in the f — B pa- 
rameter space of dipole-allowed transitions at frequencies 
/. These df/dB c^ points were identified as OWPs: Line 
narrowing and reduced sensitivity to temporal and spatial 
noise in B over a broad region of fields (closely related to 
df/dB = extrema) were found. However, to date, their 
effectiveness for reducing decoherence in the real environ- 
ment of a spin bath remains untested. 

In natural Si, 4.67% of lattice sites are occupied by the 
nuclear spin-half ^^Si isotope, rather than the spinless ^^Si. 
Flip-flopping of the ^^Si spins provides the dominant mech- 
anism of decoherence for both Si:P and SiiBi systems at low 
temperatures. The decay of the donor Hahn spin echo for 
these systems is typically fitted to exp[— t/T2 — (t/Tso)^], 
where Tsd < 1 ms characterizes the nuclear spin diffusion, 
with n c^2 — 3^ Other relaxation processes, such as those 
arising from donor-donor interactions, are represented by 
T2. Since T2 ^ Tsd, nuclear spin diffusion remains the 
main channel of decoherence at low temperatures . ^^1^^ 

In this work, we investigate the nature of the Bi-^^Si 
interaction by means of pulsed electron- nuclear double res- 
onance (ENDOR). ^^ To obtain an ENDOR spectrum, an 
EPR spin echo is detected as a function of a radio frequency 
(rf) excitation. When the rf radiation is resonant with an 
NMR transition, changes are seen in the EPR signal if the 
populations of the relevant energy levels change. Previ- 
ous ENDOR studies of Si:Bi used rf frequencies of at least 
several hundreds of MHz^i^i and thus could not probe 
the weak couplings to a spin-bath. In contrast, here, rf 
frequencies of a few MHz were used. With this approach, 
we have measured the bismuth spin bath superhyperfine 
(SHE) couplings and determined their anisotropy. 

We also present the results of cluster correlation expan- 
sion (CCE) simulations^!^ This model has been used 



with considerable success to model central spin decoher- 
ence in Si'.FA'^^^ In Ref. [2J, weak state-mixing in Si:Bi 
was investigated by simply allowing for the variation of an 
effective gyromagnetic ratio. Here we adapt the CCE simu- 
lations to include, for the first time, the strong state-mixing 
seen near the OWPs. A key finding is the demonstration of 
near-complete suppression of nuclear spin diffusion, even in 
natural Si: This occurs in extremely narrow regions, where 
TsD is in effect divergent,— in contrast to the broader ef- 
fect expected from the form of df /dB} A successful means 
of controlling decoherence is to employ isotopically en- 
riched samples^i^ i-^^1-^^1^^1^^1^-^1^^ — which can exhibit long 
T2 times up to the order of seconds ^^ Thus, the OWPs 
represent a potentially complementary technique, effective 
for both natural Si and partially enriched samples. In ad- 
dition, our work suggests that the OWPs may also be ef- 
fective in suppressing residual effects such as donor-donor 
interactions which are responsible for shortening T2 times 
from seconds to the ms timescale^ 



II. SPIN HAMILTONIAN 

We investigate a spin system with total Hamiltonian 

H = AIS ^ Hzee + ^int + ^bath- (1) 

The first term denotes the isotropic hyperfine interac- 
tion between the bismuth electronic and nuclear spins. 
For Si:Bi, the usual high-field reduction to Ising form 
AI • S ::^ AI^S^ cannot be made at the fields of inter- 
est here. The second term represents the Zeeman interac- 
tion with the external field, including the donor spins and 
summed over bath spins. 



Hzee = ^0 



(s^-SBiP-SsiY^m. 



(2) 



Here, ujq = fiB/h is the electronic Zeeman frequency 
(/i = 1.857 X 10-2^ JT-i) while fei = 2.486 x IQ-"^ and 
Ssi = 3.021 X 10~^ are the ratios of the nuclear to elec- 
tronic Zeeman frequencies for the donor and ^^Si spins, 
respectively. 

The spin-bath interaction term 



-^int — ^_^InJn^7 



(3) 



represents the SHF couplings between the donor and bath 
spins, in general of tensor form (for anisotropic couplings). 
Finally, dipolar coupling between each pair of ^^Si spins is 
represented by the bath term 



^bath — 7 ^ InD(r^rn)Im, 



(4) 



where Vnm denotes the relative position vector of bath spins 
at lattice sites n and m. Writing r^m = r for a pair of spins, 
the components of the dipolar tensor D are given by 



A,(r) 



/iO%/i' 



2 ,,2 



Sin 



r2 )' 



with /io = 47r X 10~^ NA~^ and z, j = x^y^z 



(5) 



III. ENDOR MEASUREMENTS 

The experimental ENDOR studies reported here served 
to investigate and characterize the isotropic/anisotropic 
character of the spin-bath interaction term, namely a set of 
distinct J^ values - SHF couplings - in Eq. (|3|), correspond- 
ing to occupancy of inequivalent lattice sites by ^^Si impu- 
rities. Pulsed ENDOR experiments were performed using 
the Davies ENDOR pulse sequence.^^'^^ We applied the 
pulse sequence TTmw— '?"i— TTrf— r2 — f mw— rs— TTmw— 'ra— echo, 
where the microwave (mw) frequency is chosen to excite 
one EPR transition and the rf is stochastically varied be- 
tween 2 — 12 MHz or 2 — 7 MHz to excite all nuclear spin 
transitions in this region. We used 256 ns long TTmw-pulses 
and a 128 ns long ^mw-pulse. For optimal signal-to-noise 
ratio and resolution, we used a TTrf-pulse of 10 /is. Pulse 
delays were set to n = 1 /is, r2 = 3 /is, and ra = 1.5 /is 
and a shot repetition time of 1.3 ms was employed to give 
a good signal-to-noise ratio. All experiments were carried 
out at 15 K on an E580 pulsed EPR spectrometer (Bruker 
Biospin) equipped with pulsed ENDOR accessory (E560D- 
P), a dielectric ring ENDOR resonator (EN4118X-MD4), 
a liquid helium flow cryostat (Oxford CF935), and an rf 
amplifier (ENI A-500W). We used a donor concentration 
of 3 X 10^^ cm~^ and the magnetic field was directed per- 
pendicular to the [111] plane. 

While not offering the higher frequency resolution at- 
tainable with continuous-wave ENDOR^^i^ the pulsed 
ENDOR measurements permit us to adequately constrain 
and to demonstrate the reliability of the numerical simu- 
lations. In particular, we established that isotropic cou- 
pHngs to the spin bath dominate the decoherence dynam- 
ics. While not the focus of this study, a further motivation 
is to investigate the feasibility of an alternative possibility 
for QIP: to simultaneously manipulate the ^^Si atoms as 
spin-half qubits, along with the donors^ 

Measured ENDOR spectra at / :^ 9.755 GHz are pre- 
sented in Fig. [TJ together with a list of SHF couplings. For 
the magnetic field range B c^ 0.1 — 0.6 T in Fig. (TJ there is 
significant mixing of the high-field Si:Bi energy eigenstates 
\'^s-,'^i)' The mixed eigenstates, |±,77i), correspond to 
doublets (at most) of constant m = ms + mi: 



|±,m) 



|±^,m 



.± |2 



f.±|2 



^m(^o) 



V(l^^M + 25-m2) 



T\,m±\), 

7m (wo), 



(6) 

(7) 



where ^^^(cjo) = ^^ + ^(1 + fei) and mi is an integer, 
— 5 < 771 < 5. Such mixing leads to a complex EPR spec- 
trum for bismuth with df /dB = extrema. The minima 
correspond to transitions between states corresponding to 
adjacent avoided level crossings, of which there are four. 
The disparity between the electronic Zeeman and hyper- 
fine energy scales and SHF energy scales means that the 
tensor coupling in Eq. ([3]) reduces to simpler form 



i^inM«(a^/f + A/^)5^ 



(8) 



written for coupHng to a single ^^Si at site /, where 
ai = [{a,,o,i - Ti) + 3T2cos2^^] and A = 3T^ sin ^^ cos i?^ 
with aiso,i and Ti the isotropic and anisotropic parts of the 



molecular- frame SHF tensor, respectively, and i^i the an- 
gle between the external field and the line connecting the 
bismuth site and site I. Nonsecular terms involving S^ 
and S^ can be neglected.^^ Diagonalization of the result- 
ing 2-dimensional Hamiltonian, and setting Ti = for a 
purely isotropic coupling leads to a simple expression for 
the ENDOR resonance frequency for donor level |±,77i): 



Isotropic superhyperfine couplings (MHz) 






-UJoSsi ± (^^j 7m (wo) 



(9) 



The above expression is in perfect agreement with full nu- 
merical diagonalization. The couplings in Fig. [1] were ex- 
tracted from the measured spectra by fitting to the data 
Gaussians of equal width and using Eq. ([9j) . The same ex- 
pression and a single set of couplings gave excellent agree- 
ment with data at 10 different fields. In particular, the 
observed pattern of half a dozen highest frequency ^^Si res- 
onances moving to a minimum at B c:^ 0.2 T, then increas- 
ing again, is directly attributable to mixing of the states of 
the bismuth donor: i.e., here 7m (^o) has a minimum. 

Ten out of the twelve couphngs extracted from data 
were found to be purely isotropic. The highest-field 
spectrum was measured for a range of crystal ori- 
entations and only three weak intensity lines showed 
orientation-dependent frequencies and hence anisotropy. 
Two are indicated by Xi and X2 in Fig. [TJ The cor- 
responding two couplings with nonzero anisotropy were 
found to have {aiso,Xi — 2.8, Tx^ — 2.4) MHz and 
(<^iso,X2 — 0-4, Tx2 — 2.8) MHz by fitting the more gen- 
eral form of Eq. (|9|) with nonzero T. A previous ESEEM 
(electron spin echo envelope modulation) study identified a 
single anisotropic coupling^ attributed to £^-shell (nearest 
neighbor) ^^Si. The third line we identify is fitted by cou- 
plin g co nstants consistent with the anisotropic coupling in 
Ref. [42|. For most crystal orientations, this line is masked 
by much higher intensity lines arising from isotropic cou- 
plings. 

At fields where 7m (^0) becomes small (this occurs close 
to the df /dB = minima as shown in Refs. 26.27) , Eq. (|9]) 
tends to the ^^Si Zeeman frequency (^siC^o- It is straight- 
forward to extend Eq. (|9|) to the anisotropic case and show 
that the latter statement also holds for anisotropic cou- 
plings. In effect, at these points, the donor might be said 
to approximately decouple from the bath. For example, 
for the EPR transition |12) -^ |9) (labeling the eigenstates 
|n = 1, 2, . . . 20) in increasing order of energy), 7m (^0) = 
at 5 = 157.9 mT for level |12) and B = 210.5 mT for |9). 
We note that there is however no 5-field value where both 
the upper and lower levels have 7m (^0) = 0: As we see be- 
low, this is not actually essential for complete suppression 
of spin diffusion. The actual OWP is at 5 = 188.0 mT, 
where 7-3 (cjo) = — 7_4(co'o). This is extremely close to 
where df /dB = 0, which occurs when 
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FIG. 1: (Color online) Pulsed ENDOR measured for bismuth- 
doped silicon with frequency 9.8 GHz at which ten EPR lines 
are observed, the resonance peaks due to interactions of the 
donor with ^^Si nuclei at inequivalent lattice sites. The isotropic 
superhyperfine couplings were extracted from the spectrum at 
the highest magnetic field. As the field is varied, the smooth 
lines follow the resonance positions according to Eq. (O. Solid 
and dotted lines distinguish between the two peaks observed for 
each coupling, each corresponding to one of the two donor levels 
involved in the EPR transition. Only the peaks labeled Xi and 
X2, in addition to a third pair not resolved here, were found 
to show anisotropy from performing ENDOR as a function of 
crystal orientation. 



IV. CLUSTER CORRELATION EXPANSION 
SIMULATIONS 

In order to model the full dynamics, we assume that the 
temperature and donor concentrations are low enough so 
that phonon-induced decoherence and decoherence due to 
interactions between donors are negligible. The Hahn spin 
echo decay of a central donor electron coupled to the bis- 
muth nucleus in a bath of ^^Si was calculated using the 
CCE^ Denoting the spin echo intensity by L(t)^ let Ls{t) 
be L(t) computed including only spins in some set or "clus- 
ter" of bath spins S. The quantity Ls{t) is defined as 



L{t) = \{Ls{t), 



(11) 



7_3(cJo) +7-4(^0) 



2fei 



(1 + fei) 



0. (10) where the product is over all clusters. Applying this defini- 

tion to Ls{t) and factoring out Lsit)^ an explicit form for 



the Ls{t) is obtained in terms of the Ls{t) and the Lc{t) 
in subsets C of 5, 



Coherence time, Tg„ (ms) 



Ls{t) = Lsit)/ n Lcit). 



(12) 



CdS 



The problem of calculating L{t) is reduced into indepen- 
dent components each for a distinct cluster of bath spins. 
The exact solution to L(t) is obtained if the Ls(fi) from all 
clusters are combined using Eq. (pTj) and the approxima- 
tion to L (t) up to a maximum cluster size of k is defined 



LW(t)= n Ls{t), 



(13) 



\S\<k 



which involves calculating reduced problems for all clusters 
each containing at most k spins. The L^^^'^\t) (2-cluster) 
calculation is a good approximation to L{t) when consider- 
ing only dipolar interactions in the bath affecting the spin 
echo, as these are at most a few kHz and hence perturba- 
tive compared to the SHF interactions in the MHz range 
involving the donor electron. The CCE is exact, but not 
always convergent. We calculated the 2-cluster {k = 2) ap- 
proximation to the CCE and obtained convergence for up 
to /c = 3. 

The Hahn echo sequence evolves the combined system- 
bath state to time t = 2r: 



m = 2t)) 



iHt 



"(a,®lB)e-'^n^(i = 0)), (14) 



where ax is the Pauli-X gate acting on the donor and 
1b denotes the bath identity. We assume that the time 
taken for a 7r-pulse is small compared to r. The initial 
state was written as a product of the initial donor and 
bath states, the former chosen as an equal superposition 
of states 1 12) and |9). The donor subsystem is recovered 
after tracing over the bath and the modulus of the nor- 
malized off-diagonal element of the donor reduced density 
matrix is proportional to the intensity of the echo at time 
t = 2r. The reduced problem of the 2-cluster bath was 
solved for each of the four initial 2-product bath states 
and the average intensity obtained. 2-clusters were formed 
by pairing ^^Si spins separated by up to the 3rd nearest 
neighbor distance in a diamond cubic lattice of side length 
160 A. Convergence was obtained as the lattice size and 
the separation between the two bath nuclei were extended. 
It was assumed that B was large enough to conserve the 
total ^^Si Zeeman energy. Thus, the dipolar interaction be- 
tween the two bath spins took the form of a combination 
of Ising (/f/l) and flip-flop {I^I^ + A"^^) terms. The 
Kohn-Luttinger electronic wavefunction was used to cal- 
culate the isotropic Fermi contact SHF strength with an 
ionization energy of 0.069 eV for the bismuth donor. Cal- 
culated couplings were of the same order as those obtained 
from data. The data suggest that isotropic couplings pre- 
dominate; hence anisotropic couplings were neglected and 
the simulations were largely insensitive to orientation. Fi- 
nally, we obtained the average L^^^'^\t) over 100 spatial 
configurations of ^^Si occupying 4.67% of lattice sites. The 
resulting decay curves were fitted to exp[— t/T2 — (t/Tso)^], 
obtaining T2 ^ Tsd and values of n c^ 2 — 3. 
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FIG. 2: (Color online) Suppression of Bi-^^Si spin bath de- 
coherence for the 1 12) -^ |9) EPR transition, (a) Simulated 
ENDOR and nuclear spin diffusion coherence times Tsd as a 
function of magnetic field B, showing collapse of the superhy- 
perfine couplings and a sharp increase in Tsd as the field ap- 
proaches the B = 188.0 mT optimal working point (OWP). The 
dashed line is a fit. (b) Simulated ENDOR at the B = 188.0 mT 
OWP (upper panel) and experimental spectrum at 9.755 GHz 
(lower panel), (c) Calculated donor Hahn spin echo decays from 
which coherence times in Fig. [IJa) were extracted. 



V. SUPPRESSION OF NUCLEAR SPIN 
DIFFUSION 

The results of our CCE simulations are presented in 
Fig. [21 Figure [2fa) shows the behavior around the 
B = 188.0 mT OWP, associated with the |12) -^ |9) EPR 
line. The calculated coherence time Tsd (orange dashed 
line) is superposed on a color map showing the SHF spec- 
trum: The latter shows ENDOR spectra simulated as a 
function of B^ using Eq. (J9]) and centered about the ^^Si 
nuclear Zeeman frequency. Strikingly, as B approaches the 
OWP, the "comb" of SHF lines narrows to little more than 
the width of a single line. This suggests a drastic reduc- 
tion in the value of the SHF couplings, indicating that the 
bismuth has become largely decoupled from the ^^Si spin 
bath. 

The collapse in the SHF couplings is illustrated further 
in Fig. [2jb). The lower panel shows the measured spec- 
trum at 9.755 GHz. Using our experimentally determined 
SHF couplings, the corresponding spectrum at the OWP 
is shown in the upper panel of Fig. [2fb), demonstrating 
clearly the narrowing of the spectrum [corresponding to 
the same parameters as Fig. [2fa) but at the precise field 
value of the OWP]. 

The behavior of Tsd is also quite striking and unex- 
pected: The coherence time predicted by CCE simula- 
tions increases asymptotically at the OWP by several or- 
ders of magnitude. Away from the OWP, the results agree 
well with experimentally measured values of approximately 
0.7 ms^ In Ref.|2J, in a regime of weak state-mixing, simu- 
lations using an effective gyromagnetic ratio indicated that 
Tsd was slightly reduced (by about 5%) in a regime corre- 
sponding to lower df/dB. The present study, on the other 
hand (which in contrast to Ref. 24 employed a full treat- 
ment of the quantum eigenstate mixing), shows rather an 
effect very sharply peaked about the OWP: Nuclear spin 
diffusion is predicted to be largely suppressed, but over an 
extremely narrow magnetic field range. 

Figure [2fc) shows a sample of CCE spin echo decays from 



which Tsd times were extracted, and also serves to further 
illustrate the sharp increase in Tsd- Similar suppression 
is present for other OWPs in Si:Bi. There are df/dB = 
minima for the |15) -^ |6), |14) -^ |7), |13) -^ |8), |12) -^ 
1 9), and |11) -^ |8) transitions in the frequency range 5 — 
7.5 GHz and two maxima for |12) -^ |11) and |9) -^ |8) 
close to 1 GHz. The decoupling from the spin bath is also 
expected to lead to suppression of decoherence arising from 
the interaction with a bath of donors^ 



VI. CONCLUSIONS 

In conclusion, we present measurements of the SHF cou- 
plings between a bismuth donor and a bath of ^^Si im- 
purities which suggest that isotropic couplings dominate. 
We further demonstrate the suppression of couplings at 
OWPs. Finally, the spin echo decay of the donor is calcu- 
lated as a many-body problem and sharp divergence of the 
spin diffusion time is found at an OWP. Our study moti- 
vates experimental EPR studies in the range 5 — 7.5 GHz 
corresponding to the regions of suppressed decoherence. 
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